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I 


INTRODUCTION 


The  line-of-sight  (LOS)  of  an  aberrated  optical  system  is  defined,  in 
terns  of  the  centroid  of  its  point  spread  function  (PSF).  Using  the  Fourier 
transform  relationship  between  the  PSF  and  the  optical  transfer  function 
(OTF),  it  is  also  expressed  as  the  slope  of  the  imaginary  part  of  the  OTF  at 
the  origin.  Since  the  system  OTF  is  equal  to  the  autocorrelation  of  its  pupil 
function,  the  centroid  can  also  be  written  in  terms  of  the  amplitude  and 
aberration  across  its  pupil.  Using  the  expression  in  terms  of  the  pupil 
function,  it  is  easy  to  show  that  the  LOS  obtained  by  wave  diffraction  optics 
is  identical  with  that  based  on  geometrical  optics.  It  is  shown  that,  whereas 
for  an  aberration-free  pupil,  different  amplitude  distributions  across  it  give 
the  same  LOS,  regardless  of  its  shape;  for  an  aberrated  pupil  both  the 
aberration  and  the  amplitude  distribution  affect  the  LOS.  For  a  uniform 
amplitude  distribution,  the  LOS  depends  only  upon  thg  aberration  along  the 
boundary  of  the  pupil,  i.e.,  it  is  independent  of  the  aberration  across  its 
interior. 

Next,  an  optical  system  with  aberrated  but  uniformly  illuminated  annular 
pupil  is  considered.  Its  aberration  function  is  expanded  in  terms  of  Zernike 
annular  polynomials.  It  is  shown  that  only  those  aberrations  contribute  to 
the  LOS  that  vary  with  angle  as  cos0  or  sin@,  i.e.,  various  orders  of 
coma-type  aberrations.  A  simple  expression  is  obtained  for  the  LOS  in  terms 
of  the  Zernike  aberration  coefficients.  It  is  shown  that  two  different  orders 
of  Zernike  coma  (including  tilts)  with  the  same  standard  deviation  across  the 
pupil  do  not  contribute  the  same  amount  of  LOS  error;  for  a  given  standard 
deviation,  a  higher-order  Zernike  coma  gives  a  larger  LOS  error.  Finally, 
some  numerical  examples  are  considered.  In  particular,  the  PSF's  aberrated  by 
primary  and  secondary  classical  coma  are  discussed.  It  is  shown  that  such 
aberrations  shift  the  peak  and  centroid  of  the  PSF  (which  are  coincident  in 
the  aberration- free  case)  by  significantly  different  amounts.  Similar  results 
are  obtained  for  annular  pupils  with  radially  symmetric  illumination.  These 
results  are  illustrated  by  considering  Gaussian  illumination  of  an  aberrated 
annular  pupil.  Numerical  results  for  both  circular  and  annular  pupils  are 
also  given.  In  particular,  it  is  shown  that  two  PSF's  aberrated  by  the  same 


amount  of  primary  and  secondary  coma  have  identical  centroids  in  the  case  of 
uniform  illumination,  but  different  centroids  for  Gaussian  illumination. 

The  results  obtained  here  are  applicable  to  both  imaging  systems  as  well 
as  laser  transmitters.  For  a  certain  point  object,  the  amplitude  distribution 
across  the  pupil  of  an  imaging  system  will  generally  be  uniform.  However, 
across  a  laser  beam,  the  amplitude  distribution  is  often  given  by  a  Gaussian 
function.  Moreover,  laser  beams  are  often  polarized.  If  the  polarization  is 
linear,  then  the  scalar  treatment  of  diffraction  considered  in  this  paper 
leads  to  a  linearly  polarized  diffracted  wave.  If  the  beam  is  circularly  or 
elliptically  polarized,  then  the  diffraction  of  the  two  orthogonal  polariza¬ 
tion  components  may  be  treated  separately.  The  total  diffracted  field  is 
obtained  by  taking  a  vector  sum  of  the  two  diffracted  components. 

II.  THEORY 

The  PSF  of  an  optical  imaging  system,  i.e.,  the  irradiance  distribution 
of  the  image  of  a  certain  point  object,  according  to  wave  diffraction  optics 
is  given  by1 

I(x,y)  »  (1/A2R2)|  //  P(u,v)  exp(-2iri(xu  +  yv)/AR]  du  dv|2,  (la) 

S 

where  X  is  the  wavelength  of  the  object  radiation  and  i  »  /-l.  P(u,v)  is 
the  pupil  function  of  the  system  corresponding  to  this  point  object,  and  if 
A(u,v)  and  W(u,v)  represent  the  amplitude  and  aberration  of  the  light  wave  at 
a  point  (u,v)  on  its  exit  pupil,  then 

P(u,v)  *  A(u,v)  expt27TiW(u,v)/X] ,  inside  the  pupil, 

*  0,  outside  the  pupil.  (lb) 

The  integration  in  Eq.  (la)  is  carried  out  over  the  clear  or  illuminated 
region  of  the  pupil.  S  represents  the  area  of  this  region. 

The  aberration  function  W(u,v)  represents  the  deviation  of  the  optical 
wavefront  at  the  exit  pupil  from  a  spherical  wavefront,  called  the  reference 
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sphere,  measured  along  a  ray  passing  through  the  point  (u,v).  The  aberration 
W(u,v)  is  considered  positive  if  a  ray  from  the  point  object  passing  through 
the  point  (u,v)  has  to  travel  a  Longer  optical  path  in  reaching  the  reference 
sphere  than  a  reference  ray  which  passes  through  the  center  of  the  exit 
pupil.  The  optical  wavefront  and  the  reference  sphere  pass  through  the  center 
of  the  exit  pupil.  The  reference  sphere  has  a  radius  of  curvature  R.  Its 
center  of  curvature  is  generally  chosen  to  be  at  the  Gaussian  image  of  the 
point  object.  The  center  of  curvature  may  also  be  chosen  such  that  the 
variance  of  W(u,v)  across  the  pupil  is  minimized.  In  either  case,  the  center 
of  curvature  defines  the  origin  of  the  (x,y)  image  plane,  which  is  parallel  to 
the  (u,v)  plane.  The  line  joining  the  centers  of  the  exit  pupil  and  the 
reference  sphere  defines  a  reference  axis.  In  a  system  with  an  axis  of  rota¬ 
tional  symmetry,  it  may  be  considered  as  the  reference  axis. 

The  LOS  for  the  point  object  as  perceived  by  the  optical  system  is 
determined  by  the  centroid  of  its  PSF.  If  <x>  and  <y>  represent  the  coordi¬ 
nates  of  this  centroid,  then 


<x>  ■  E  1  //  x  I(x,y)  dx  dy, 


<y>  *  E  //  y  I(x,y)  dx  dy, 


where 


E  ■  SS  I(x,y)  dx  dy, 

—GO 

is  the  total  power  (energy)  in  the  image.  Applying  Parseval's  theorem  to 
Ea.  (1),  we  find  that  E  is  also  given  by 


E  m  IS  |P(u,v) | 2  du  dv 
S 

■  //  l(u,v)  du  dv  , 

S 


where 


l(u,v) 


A2(u,v) 


(4) 


is  Che  irradiaace  at  a  pupil  point  (u,v).  Equations  (3a)  and  (3b)  represent 
conservation  of  energy;  i.e.,  the  total  energy  in  the  image  is  equal  to  the 
total  energy  in  the  pupil.  If  <a>  and  <$>  represent  the  angular  LOS,  they 
are  given  by 

<a>  3  <x>/R  (5a) 


and 

<g>  *  <y>/R. 


(5b) 


where  we  have  assumed  that  the  angles  are  small  so  that  they  are  approximately 
equal  to  their  tangents. 


The  OTF  of  an  imaging  system  is  equal  to  the  Fourier  transform  of  its 
PSF.  Thus  if  t(£,  ri)  represents  the  OTF  corresponding  to  a  spatial  frequency 
(5,T}),  then 


00 

T(S.n)  *  E_l  //  l(x,y)exp[2iri(£x  +  ny)]  dx  dy. 

-CD 


(6) 


Differentiating  both  sides  of  Eq.  (6)  with  respect  to  £  and  evaluating  the 
result  at  £  3n  ”  0,  we  find  Chat 


Similarly, 


<x> 


(7a) 


<y> 


l 

2ttT 


(7b) 


Thus  the  centroid  of  the  PSF  of  an  optical  system  is  given  by  the  slope  of  its 
corresponding  OTF  at  the  origin.  However,  since  <x>  and  <y>  are  real, 


only  Che  slope  of  the  imaginary  part  of  the  OTF  at  the  origin  contributes  to 
Che  centroid.2  Thus,  we  may  write 


and 


(8a) 


(8b) 


The  OTF  is  also  given  by  the  autocorrelation  of  the  pupil  function  as 
may  be  seen  by  substituting  Eq.  (1)  into  Eq.  (6). 

Thus, 


T(£,n)  -  E-1  If  P(u,v)  P*(u  -  AR£,v  -  ARn)  du  dv,  (9) 

l 

where  *  indicates  a  complex  conjugate,  and  I  is  the  region  of  overlap  of 

two  pupils  centered  at  (0,0)  and  (AR£,ARrj).  Substituting  Eq.  (9)  into  Eq.  (8), 

we  obtain3 

<x>  =■  -(AR/2TTE)  IS  Imjp(u,v)-^~y,v)j  du  dv  ,  (10) 

and  a  similar  equation  for  <y>.  Substituting  Eq.  (lb)  into  Eq.  (10),  we  obtain 

<x>  *  //  I(u,v)  ~  ^U1V^  du  dv  .  (11a) 

E  s  du 

Similarly, 

<y>  *  ~  If  I(u,v)  ^  du  dv  .  (lib) 

Since  R(3W/3u)  and  R(3W/3v)  represent  the  ray  aberrations?  i.e.,  the  image- 
plane  coordinates  of  a  ray  passing  through  the  pupil  point  (u,v),  Eq.  (11) 
shows  that  the  centroid  of  the  PSF  according  to  wave  diffraction  optics  is 
identical  with  that  according  to  ray  geometrical  optics.5  (The  PSF  given  by 
ray  geometrical  optics  is  called  the  spot  diagram). 


-9- 


From  Eq.  (11),  we  also  note  Chat  amplitude  variations  across  the  pupil 
affect  the  LOS  only  if  it  is  aberrated.  In  the  absence  of  aberrations,  the 
PSF  centroid  lies  at  the  center  of  curvature  of  the  reference  sphere  regard¬ 
less  of  the  shape  of  the  pupil  and  the  amplitude  distribution  across  it.  This 
may  also  be  seen  from  Eqs.  (1)  and  (3).  From  Eq.  (1),  we  note  that  if 
W(u,v)  ■  0,  then  l(x,y)  ■  l(-x,-y).  Hence  the  symmetry  of  the  aberration- 
free  PSF  yields  its  centroid  at  the  origin  of  the  (x,y)  image  plane.  Similar¬ 
ly,  since  the  aberration-free  OTF  is  real  (see  Eq.  (9)),  Eq.  (8)  also  gives 
the  centroid  at  the  origin. 

It  should  be  noted  that  the  peak  value  of  an  aberrated  PSF  may  or  may 
not  lie  at  its  origin,  depending  on  the  magnitude  and  the  type  of  aberration, 
whether  or  not  the  amplitude  across  the  pupil  is  uniform.  However,  tne  peak 
value  of  an  aberration-free  PSF  always  lies  at  its  origin,  regardless  of  the 
amplitude  variation  across  the  pupil.  This  may  be  seen  from  Eq.  (1).  If  we  let 

f(u,v)  »  P(u, v )  exp[-2lti(xu  +  yv)/AR],  (12a) 


then,  using  Holder's  inequality,  we  find  that 


IS 


. 


i //  f(u,v)  du  dv| 2 

I(x,y)  .  _ S _ 

I  0,0  //  | f (u,v) | 2du  dv 

S 

<  1  . 


(12b) 


m 


Thus,  both  the  centroid  and  the  peak  value  of  an  aberration-free  PSF  lie  at 
its  origin,  regardless  of  the  amplitude  variations  across  the  pupil. 

Equations  (3),  (8),  and  (11)  give  the  LOS  of  the  system  in  terms  of  its 
PSF,  OTF,  and  the  aberration  function.  In  practice,  given  an  imaging  system, 
the  most  convenient  expression  to  use  would  be  Eq.  (3),  since  the  PSF  can  be 
measured  easily  by  using  a  photodetector  array.  In  optical  design  and 
analysis,  the  simplest  way  to  obtain  the  LOS  would  be  to  use  Eq.  (11),  since 
the  aberrations  must  be  calculated  even  if  the  other  two  expressions  were 
used.  Thus  one  may  trace  rays  all  the  way  up  to  the  image  plane  and  determine 


I--  V.-.V.A 
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the  centroid  of  the  ray  distribution  in  this  plane  with  appropriate  weighting 
I(u,v)  of  each  ray. 


If  the  illumination  across  the  pupil  is  uniform,  e.g.,  if 


A(u,v)  =  Aq  , 


and 


l(u,v) 


and,  therefore. 


then  Eq.  ( 1 1 )  reduces  to 


h  • 


SIo  - 


and 


_  _  R  rr  3W(u,v)  ,  . 
<x>  *■  —  /  J  — 5 — * —  du  dv 

s  s  9u 


„  „  R  rr  3W(u,v)  , 

<y>  s  !J  ~  8v  du  dv  • 

u 


(13a) 


(13b) 


(13c) 


(14a) 


(14b) 


Using  Stokes  theorem,  the  surface  integral  in  Eq.  (14)  involving  the 
derivative  of  the  aberration  function  can  be  written  in  terms  of  its  line 
integral  along  the  curve  bounding  the  surface.  Thus,  we  may  write6 


<x>  =  (R/S)  f  W(u,v)  u  •  ds  (15a) 

and 

<y>  *  (R/S)  f  W(u,v)  v  •  ds  ,  (15b) 

A  A 

where  u  and  v  are  unit  vectors  along  the  u  and  v  axes,  respectively,  ds 
represents  an  element  of  arc  length  vector  of  the  curve  bounding  the  pupil. 

It  is  evident  from  Eq.  (15)  that  in  the  case  of  an  aberrated  but  uniformly 
illuminated  pupil,  the  centroid  of  the  PSF  can  be  obtained  from  the  value  of 
the  aberration  function  only  along  the  perimeter  of  the  pupil.  Accordingly, 
in  that  case,  to  calculate  the  centroid  the  knowledge  of  the  aberration  across 
the  interior  of  the  pupil  is  not  needed. 
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III.  APPLICATION  TO  SYSTEMS  WITH  ANNULAR  PUPILS 
A.  Uniform  Illumination 

Consider  an  imaging  system  with  an  aberrated  but  uniformly  illuminated 
annular  pupil  of  inner  and  outer  radii  of  ca  and  a,  respectively,  where 
0  <  e  <  1.  The  area  of  the  pupil  is  given  by 

S  *  ir(l-e2)a2  .  (16) 

Owing  to  the  circular  boundary  of  the  pupil,  it  is  convenient  to  use  plane 
polar  coordinates  (h,6),  where 


and 


so  that 


and 


u  ■  h  cos© 
v  ■  h  sind, 


0  ■  tan  l(v/u) 


h  -  (u2  ♦  v2) 


1/2 


(17a) 

(17b) 


(17c) 
( 17d) 


Moreover,  0  0  <  2n  and  ea  _<  n  _<  a.  A  vector  element  of  a  circular  arc  with 

a  center  of  curvature  at >(0,0)  and  passing  through  a  point  (u,v)  is  given  by 

ds  m  (u  d0,  v  d0),  (18) 

Let  the  aberration  function  in  polar  coordinates  be  W(h,0;e).  Substituting 
Eq.  (18)  into  Eq.  (15),  we  obtain 


<x>  -  [R/ir(l-e2)aJ 


2» 

/  [w(a,0;e)  -  e  W(ea,0;e)]  cos0  d0 
0 


(19a) 


and 


2TT 

<y>  *  [R/it( l-e2)a]  /  [W(a,9;e)  -  e  W(ea,0;e)]  sin0  d0 

0 


(19b) 


LeC  us  expand  Che  aberration  function  W(h,0;e)  in  terms  of  a  complete 

:  Zernixe  annular  polynomials  Rm(p;e) 

n 

orthogonal  over  the  annular  pupil7,  where 


set  of  Zernixe  annular  polynomials  R  (p;e)  cosm0  and  R  (p;e)  sinm0  which  are 

n  n 


p  *  h/a. 


(20) 


Thus,  we  may  write 


W(b,9;e)  *  I  J  e  /2(n+l)  R*(p;e)  (c  cosmO  +  s  sinm0), 

n  n  m  a  nm  am 

n*0  m**0 

where  n  and  m  are  positive  integers  (including  zero),  n-m  >_  0  and  even. 


(21) 


e  -  1//T,  ■  -  0 


•  1 


i  0, 


(22) 


and  c  and  s  are  the  expansion  coefficients.  Rote  that 
nm  nm 


nO 


0. 


The  radial  polynomials  obey  the  orthogonality  relation 


/  R>;c)  R*,(o;e)  P  dp  -  , 


(23) 


(24) 


where  6  ,  is  a  KronecXer  delta.  The  magnitude  of  each  expansion  coefficient 

nn 

(except  Cqq)  represents  the  standard  deviation  of  the  corresponding 
aberration  term  across  the  annular  pupil. 
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Substituting  Eq.  (21)  into  Eq.  (19),  and  noting  the  orthogonality  of 
trigonometric  functions,  we  obtain 


.00 

<X>  =  [R/(l-e2)a]  l'  / 2 ( n+ 1 )  [R1(l;e)  -  e  Rl(e;e)J  c  .  (25a) 

n-1  Q  n  nl 

and 

00 ,  _ 

<y>  3  [R/ ( l-e2)a]  £  /2(n+l)  {R^lje)  -  e  Rl(e*,e)J  s  ,  (25b) 

,  n  n  ni 

n*l 

where  a  prime  on  the  summation  sign  indicates  a  summation  over  odd  integral 

values  of  n.  Thus  the  only  aberrations  that  contribute  to  the  LOS  error 

are  chose  with  mol.  Aberrations  of  the  type  R^(p;e)  cos0  contribute  to 

1  a 

<x>  and  those  of  the  type  &  (p;e)  sin0  contribute  to  <y>,  This  is  also 

a  j_ 

evident  from  the  symmetry  of  the  aberrations.  Since  R  (p;e)  consists  of  terms 

n  n  2  1  Q 

in  p  ,  p  , and  p,  therefore,  for  example,  R*(p;e)cos0  is  symmetric  in 

v  but  not  in  u.  Hence  the  PSP  is  symmetric  in  y  as  may  be  seen  from  Eq.w (1) 

noting  that  che  amplitude  across  the  pupil  is  uniform..  Accordingly,  <y>  ■  0 

for  this  aberration. 

It  is  well  known  that,  for  stall  aberrations,  the  Strehl  ratio  for  an 
optical  system  depends  on  the  variance  of  its  aberration*.  Therefore,  two 
aberration  terms  with  m  ■  l,  but  different  values  of  n,  affect  the  Strehl 
racio..in  the  same  way  if  their  coefficients  are  equal  in  magnitude.  However, 
it  is  evident  from  Eq.  (25)  that  their  contribution  to  the  LOS  error  is  not 
the  same;  an  aberration  of  higher  order  contributes  a  larger  error  for  the 
same  value  of  the  coefficient.  Hence,  in  tolerancing  an  optical  system,  one 
must  be  careful  in  allocating  equal  standard  deviation  to  two  aberration  terms 
that  also  contribute  to  the  LOS  error. 

In  the  case  of  a  circular  aperture  (e  ■  0), 

R1(1;0)  -  1  .  (26) 

n 
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Therefore,  Eq.  (25)  simplifies  to 


<x>  »  (R/a)  y'  /2(n+l)  c 

nx 

n*  l 

and 

00  f 

<y>  -  (R/a)  V  /2(n+l)  s  . 

n-1  nl 

We  aote  that  as  in  the  case  of  an  annular  pupil,  two  aberration  terms  with 
equal  standard  deviation  but  different  order  n  do  not  contribute  equally  to 
the  LOS  error.  For  a  given  standard  deviation,  a  higher-order  aberration 
gives  a  larger  LOS  error  compared  to  a  lower-order  aberration. 

B.  Radially  Symmetric  Illumination 

Let  A(h)  and  1(h)  describe  the  radially  symmetric  amplitude  and 
irradiance  distributions  across  the  aberrated  annular  pupil,  where 

1(h)  -  A*(h).  (28) 


(27a) 


(27b) 


In  polar  coordinates,  Eq.  (11)  for  the  LOS' can  be  written 


and 


a 

/ 

Ca 


a 

/ 

Ca 


2ir 

/ 

0 


2ir 

/ 

0 


1(h)  cos0 


1(h)  sind 


3K(.fr|9;e)  _  sine  3wOhiL£l1  h  dh  d0 

dh  h  36  J  a  aD 


3W(h,0;c)  .  cos6  3w(h,0;e)l  .  lA 

- aS —  *  — - se^—J  h  dh  de 


where 

a 

E  *  2tt  /  1(h)  h  dh  d0  . 

Ca 


(29a) 


(29b) 


(30) 
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We  now  expand  Che  aberration  function  in  terms  of  anhular  polynomials 
S^(p;e)  cosra0  and  S^(p*,e)  sinm0  that  are  orthogonal  over  the  radially 
symmetric  illuminated  annular  pupil7.  Thus  we  write 


a* 

W(h,0;  e)  “I  F  e  /2(n+l)  Sm(p;e)(c  cosm0  +  s  ainm0),  (31) 
__  _  m  n  nm  nm  ■ 


n*0  mM0 


where 


s“(p;e) 


(n-m)/2 

R®(p;e)  -  l  (n-2i+l)<R®(p; e)  s;_,.(p;e)>  s”_,.(p;e) 
i>i 


n-2i 


(32) 


is  a  radial  polynomial  with  properties  similar  to  those  of  R™ (p;e).  It  should 

n 

be  evident  that  the  aberration  coefficients  c  and  a  in  Eq.  (31)  are 

nm  nm 

different  from  those  in  Eq.  (21).  In  Eq.  (32) 


<R“(p;e)  S^_2i(p;e)> 


12  it  1  2tt 

/  /  R*(p;e)  s“  ,.(p;e>  A(.p)p  dp/  S  S  A(.p)  p  dp, 
p  0  C  o  (33) 


and  M™  is  a  normalization  constant  such  that 
n 


1  1  l 

/  S®(p;e>  S®,(p;c)  A(p)  p  dp  /  /  A(p)  p  dp  -  —  ^  6qq, 


(34) 


Substituting  Eq.  (31)  into  Eq.  (29)  we  find  that 


<x> 


TTaR 

E 


°°l  1  ^ 

l'  /2(n+l)  cq1  /  l(ap)  tp  S*(p;c) J  dp 


n«l 


(35a) 


and 


00|  i  n  , 

<y>  -  I  /2(n+l)  snL  /  l(ap)  |-  (p  s*(p;e)]  dp.  (35b) 


n*l 


16- 


Once  again  Che  only  aberraciona  that  contribute  to  the  LOS  error  are 

thoae  with  m  -  l.  Aberrations  of  the  type  S^(p;e)  coa0  contribute  to  <x>  and 

1  '  n 

thoae  of  the  type  Sn(p;e)  sin8  contribute  to  <y>,  aud  this  can  be  explained 
from  the  symmetry  properties  of  the  aberrated  PSF  in  a  manner  similar  to  the 
discussion  following  Eq.  (25). 

IV.  NUMERICAL  EXAMPLES 

A.  Uniform  Illumination 


Using  polar  coordinates 


and 


x  ■  r  cos<|> 
y  ■  r  sin<{>, 


(36a) 

(36b) 


it  can  be  shown  that  for  an  aberration-free  optical  system  with  a  uniformly 
illuminated  annular  pupil,  Eq.  (1)  reduces  to 


I(r;e) 


(1-e2)2 


2J,  (trr  )  2J.  (irer  )" 

— i - 2_  -  e2  — - - — 


Trr 


irer 


l(0;e) 


Here- 


2  _i_  ..2x1/2 


r  ■  (xz  +  y  ) 


(37) 


(38) 


is  the  radial  distance  of  a  point  (x,y)  in  the  image  plane  from  its  origin, 
and 

rg  -  r/XF  (39) 

is  a  scaled  radial  distance  of  an  image  point.  J^(*)  is  the  first-order 
Bessel  function  of  the  first  kind,  and 


F  -  R/2a 


(40) 
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ia  the  f-number  of  the  system.  The  central  value  of  the  PSF  ia  given  by 

l(0;e)  -  ES/A2R2,  (41) 

where  the  power  E  in  the  image  is  given  by  Eq.  (13c). 

The  line  joining  the  center  of  the  pupil  and  the  centroid  of  the  PSF  for 
a  given  point  object  defines  the  LOS  of  the  optical  system  in  its  image  space 
for  that  point  object.  Since  the  aberration-free  PSF  is  circularly  symmetric 
about  the  origin  of  the  (x,y)  plane,  its  centroid  also  lies  at  this  origin, 
i.e.,  it  lies  at  the  center  of  curvature  of  the  reference  sphere  with  respect 
to  which  the  aberration  for  the  point  object  under  consideration  is  zero. 
Hence,  the  line  joining  the  origins  of  the  (u,v)  and  (x,y)  planes  defines  the 
LOS. 


When  aberrations  are  introduced  into  the  system  such  that  the  centroid 

of  its  PSF  shifts,  the  pos.ition  of  the  point  object  as  perceived  by  the  system 

changes,  i.e.,  there  is  a  LOS  error.  We  have  shown  that  in  the  case  of 

optical  systems  with  pupils  having  circular  boundaries  and  uniform  or  radially 

symmetric  amplitude  distributions,  the  only  aberrations  that  contribute  to  the 

LOS  error  are  of  the  type  R1(p;E)cos0  and  R*(p;e)sin0.  Since  the  radial 
1  **  n  Q  q_  2 

polynomials  RQ(p;e)  consist  of  terms  in  p  ,  p  , ...,  and  p,  with  their 
coefficients  varying  with  e,  there  is  no  loss  of  generality  if  we  consider 
aberrations  of  the  type  pacos0  (or  pnsin0),  where  n  is  an  odd  integer,  to 
determine  their  contribution  to  the  LOS  error.  Thus,  we  consider  an  aberration 


W(h,0)  ■  Wq  (h/a)n  cos0,  ea  £  h  £  a 


(42a) 


■  W  p°  cos0,  e  <  p  <  t,  (42b) 

n  —  — 

where  is  the  peak  aberration  at  the  edge  of  the  pupil  relative  to  a  value 
of  zero  at  its  center.  Since  <y>  and  correspondingly  <$>  are  both  zero 
for  this  type  of  an  aberration,  we  shall  not  explicitly  so  state  from  now  on. 

We  shall,  for  example,  refer  to  the  centroid  as  simply  <x>. 
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Substituting  Eq.  (42)  into  Eq.  (20),  we  obtain 


<x> 


2W  P 
n 


(n-l)/2 

l 

i-0 


The  corresponding  angular  LOS  error  is  given  by 


(43) 


where 


<o> 


2(W  /D) 
n 


u-o/2 

I 

i*0 


(44) 


D  ■  2a 


(45) 


is  the  outer  diameter  of  the  annular  pupil.  We  note  from  Eq.  (44),  for 

example,  that  for' circular  pupils  (c  ■  0),  one  wave  (W  *  IX)  of  aberration  of 

n 

the  type  given  by  Eq.  (42)  produces  an  angular  LOS  error  of  2X/D.  This  is 
quite  large  considering  that  the  angular  radius  of  the  (aberration-free)  Airy 
disc-is  only  1.22  X/D. 

It  is  interesting  to  note  that  when  e  *  0,  the  LOS  error  depends  only 

on  the  value  of  W^  but  not  on  n,  the  power  of  p  in  Eq.  (42b).  This  is 

consistent  with  our  earlier  observation  (following  Eq.  (15)]  that,  for  a 

uniformly  illuminated  pupil,  the  centroid  of  its  aberrated  PSF  depends  only  on 

the  aberration  along  its  perimeter.  In  the  case  of  a  circular  pupil,  the. 

variation  of  the  aberration  along  its  perimeter  as  given  by  Eq.  (42)  is  the 

same  for  different  orders  of  the  aberration.  Hence,  for  a  given  value  of 

W  ,  even  though  aberrations  such  as  tilt  (n  ■  1),  primary  coma  (n  •  3), 
n 

secondary  coma  (n  ■  5),  etc.,  which  are  completely  different  from  each  other 
across  the  interior  of  the  pupil  and,  therefore,  give  completely  different 
PSF* s,  nevertheless  give  the  same  centroid  since  the  aberrations  are  identical 
on  the  perimeter  of  the  pupil.  This  is  not  true  for  an  annular  pupil,  in 
which  case,  although  the  aberration  along  the  outer  perimeter  is  the  same  for 
different  values  of  n,  it  is  different  along  the  inner  perimeter.  Hence,  for 
a  given  value  of  Wq,  an  annular  pupil  gives  aberrated  PSP's  with  different 
centroids  for  different  orders  of  the  aberration. 
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We  may  add  Chat  similar  observations  hold  for  balanced  aberrations 

represented  by  Zernitce  polynomials  R*(p;e)  cos0  and  R^(p:e)sin6.  For 

a  .  n 

example,  for  a  given  value  of  W  ,  aberrations  W  RL(p;0)cos6  across  a 

n  n  n  | 

circular  pupil  give  PSF's  with  the  same  centroid  since  R*(1;0)  ■  1,  i.e. , 
since  the  aberration  along  the  perimeter  of  the  pupil  is  independent  of  the 
aberration  order  n.  In  the  case  of  an  annular  pupil,  R^(l;e)  as  well  as 
R^(e;e)  depend  on  the  value  of  n.  Hence,  PSF's  with  different  centroids  are 
obtained  for  different  aberration  orders  n. 

Wavefront  Tilt 

If  we  let  n  *  1  in  Eq.  (42),  the  aberration  is  simply  a  tilt  of  the 
optical  wavefront  through  the  center  of  the  pupil.  The  PSF  shifts  and  its 
centroid  moves  from  (0,0)  to 

<x>  -  2WlF  (4oa) 

corresponding  to 

<0t>  -  2Wl/D.  (46b) 

y 

Thus,  for  example,  one  wave  of  wavefrcnt  tilt,  produces  an  angular  LOS  error 
of  2\/0,  regardless  of  the  value  of  e.  Of  course,  in  this  case,  the 
position  of  the  peak.  value  of  the  PSF  is  coincident  with  the  position  of  its 
centroid. 

Primary  Coma 

If  we  let  n  ”  3,  in  Eq.  (42)  the  aberration  obtained  is  called  classical 
primary  coma.  The  centroid  in  this  case  is  given  by 

<x>  -  2W3P(1  +  e2).  (^7) 

Thus,  for  a  given  value  of  W^,  the  centroid  for  an  annular  pupil  shifts  by  a 
factor  of  (1  +  e2)  larger  than  that  for  a  circular  pupil. 
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For  small  values  o£  W^,  Che  peak  value  of  Che  aberraced  PSF  occurs  aC 
a  point  such  Chat  if  Che  aberration  is  measured  with  respect  Co  a  reference 
sphere  centered  at  this  point,  the  variance  of  the  aberration  across  Che 
annular  pupil  is  minimum.  From  the  properties  of  the  Zernilce  annular  polyno¬ 
mials,  we  find  that  the  polynomial  R„(p;e)cos0  gives  the  optimum  combination 
of  p  cos6  and  pcoad  terms  leading  to  a  minimum  variance.  Since 

a‘(Pie)  -  3(Uca)p’  -  UHsOtsSz-  ,  (48) 

3  (l-e2)[(l+e2)(l+4c2+e1*)] 1 

we  note  that,  for  small  values  of  W^,  the  peak  value  of  the  aberrated  PSP 
occurs  at 

X  -  4W_F  (1  +  e2  +  e*)/3(l  +  e2),  (49) 

m  3 

where  the  subscript  m  refers  to  the  point  corresponding  to  minimum  aberration 
variance.  From  the  form  of  the  aberration,  it  is  understood  that  y^  -  0. 

Thus,  an  amount  of  primary  u-coma  shifts  the  centroid  and  peak  of  the  PSP 
by  different  amounts,  the  movement  of  the  peak  being  2(1  ♦  e2  +  e%)/3(l  +  C2)2 
of  the  movement  of  the  centroid.  As  an  example,  a  circular  pupil  aberrated  by 
a  quarter  wave  of  primary  u-coma  (W^  ■  X/4)  gives  an  aberrated  PSP  with  a 
centroid  at  <x>  *  XP/2  and  a  peak  value  at  x„.  •  XF/3. 

m 

For  large  values  of  W3,  the  peak  of  the  aberrated  PSP  does  not  occur  at 

the  point  corresponding  to  minimum  aberration  variance*®.  For  example,  in  the 

case  of  circular  pupils,  the  peak  lies  approximately  at  the  point  corresponding 

to  W,  -  (2/3)W_  only  when  W,  <  0.7X.  For  larger  values  of  W-,  the  peak  occurs 

closer  to  the  origin  than  the  point  corresponding  to  minimum  aberration  variance. 

For  W  >  1.6X,  the  distance  of  the  peak  from  the  origin  does  not  increase  mono- 
3  - 

tonically,  but  fluctuates  as  W3  increases.11  Since,  according  to  Eq.  (47),  the 
distance  of  the  centroid  increases  linearly  with  W3,  it  is  clear  that  the  separa¬ 
tion  between  the  locations  of  the  centroid  and  peak  increases  as  W3  increases. 
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If  we  let  4>  -  0  and  A(h)  ■  Aq  in  Eq,  (A6),  we  find  that,  along  the  x-axis, 
the  PSF  aberrated  by  primary  coma  of  the  type  W3p  cos0  can  be  written 

1 

I(x;e)  ».[I(0;e)/(l-e2)2)  [  /  Jrt(wB)  dt]2,  (50) 

e2  0 

where 

B  -  (2tW,  -  x 

J  8 

x  ■  x/AF, 

8 

and  is  in  units  of  X.  Figure  1  shows  how  l(x;e)  normalized  by  the 

aberration-free  central  value  l(0;e)  given  by  Eq.  (41)  varies  with  x  for  several 

typical  values  of  varying  from  0  to  2\,  and  e  ■  0  and  e2  ■  0.5.  Figure  2 

shows  how  the  irradiance  I  (W~;e)  at  x  ,  the  peak  irradiance  I  (W~;e)  and  the 

mo  m  po  .. 

irradiance  I  (W,;e)  at  <x>  vary  with  W,.  Figure  3  shows  how  x  ,  x  (the  point 

C  j  j  ID  p 

at  which  peak  irradiance  occurs),  and  <x>  vary  with  Wy  The  observations  made 

above  about  the  PSF's  aberrated  by  primary  coma  are  evident  from  these  figures. 

Several  typical  values  of  x  ,  x  ,  and  <x>  and  the  corresponding  irradiances 

m  p 

1,1,  and  I  are  noted  in  Table  1,  where  the  numbers  without  parentheses  are 
m  p  c  r 

for  e  m  0  and  those  with  parenthesis  are  for  e2  ■  0.5.  The  aberrated  central 
irradiance  1(0)  is  also  given  in  this  Table.  The  irradiance  values  1(0)  and 
1^  are -the  Strehl  ratios  calculated  for  primary  and  balanced  primary  coma, 
respectively,  in  an  earlier  paper.10 

Secondary  Coma 

If  we  let  n  ■  5  in  Eq.  (42),  the  aberration  obtained  is  called  classical 
secondary  coma.  The  centroid  in  this  case  is  given  by 

<x>  -  2W5F(1  +  e2  +  e1* ) .  (53) 


(51) 

(52) 
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Figure  1.  PJF  I(x  ; c)  for  several  typical  values  of  primary  coma  abarra- 

9 

tion  in  units  of  X.  The  amplitude  A(u,v)  across  the  pupil  is 
uniform.  The  PSF's  are  normalized  by  the  aberration-free  central 
value  given  by  Eq.  (41).  represents  x  in  units  of  XF. 

(a)  e  -  0,  (b)  e2  -  0.5.  * 
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Table  1.  Typical  values  of  x  ,  x  ,  and  <x>  in  units  of  XF,  and  the 

ra  p 

corresponding  irradiances  1,1,  and  I  in  units  of  the  aberration-free 

n  p  c  j 

central  irradiance  for  PSF's  aberrated  by  primary  coma,  W(h,6)  ■  W,p  cosO. 
The  units  of  are  X.  The  aberrated  central  value  1(0)  is  also  given  here. 
The  amplitude  A(u,v)  across  the  pupil  is  uniform.  The  numbers  without 
parentheses  are  for  a  circular  pupil  (e  *  0)  and  those  with  parentheses are 
for  an  annular  pupil  with  £  ■  0.5. 


u 

3 

X 

m 

X 

P 

<x> 

I 

m 

I 

P 

I 

c 

1(0) 

0 

0 

0 

0 

1.0000 

1.0000 

1.0000 

1.0000 

(0) 

(0) 

(0) 

(1.0000) 

(1.0000) 

(1.0000) 

(1.0000) 

0.5 

0.67 

0.66 

1.00 

0.8712 

0.8712 

0.6535 

0.3175 

(0.78) 

(0.78) 

(1.50) 

(0.9283) 

(0.9283) 

(0.0524) 

(0.0403) 

1.0 

1.33 

1.30 

2.00 

0.5708 

0.5717 

0.1445 

0.0791 

(1.56) 

(1.55) 

(3.00) 

(0.7410) 

(0.7412) 

(0.1357) 

(0.0319) 

1.5 

2.00 

1.80 

3.00 

0.2715 

0.2844 

0.0004 

0.0618 

(2.33) 

(2.32) 

(4.50) 

(0.5063) 

(0.5064) 

(0.0139) 

(0.0010) 

2.0 

2.67 

1.57 

4.00 

0.0864 

0.1978 

0.0061 

0.0341 

(3.117 

(3.07) 

(6.00) 

(0.2936) 

(0.2946) 

(0.0160) 

(0.0000) 

2.0 


The  variance  of  an  aberration  of  the  type  p^cosd  is  reduced  if  an  amount 
-t(l+e2+e'*+e6)/2(l+e2)]W5  of  pcosd  aberration  is  introduced. 

Hence,  the  corresponding  value  .of  x  is  given  by 

xm  »  W,F(l+e2+e*+e6)/(l+e2).  (54) 

Thus,  for  example,  in  the  case  of  a  circular  pupil,  the  variance  is  reduced 

by  a  factor  of  4  if  ■  -0.5W Accordingly,  for  small  values  of  W^,  the 

peak  of  the  corresponding  aberrated  PSF  occurs  at  x  *  WCF,  while  the  centroid 

m  j 

occurs  at  <x>  *  2W^F.  The  aberrated  PSF  along  the  x-axis  in  this  case  is 
given  by  Eq.  (50)  where 

B  -  (2t2  W5  -  xg)t1/2,  (55) 

and  is  in  units  of  X.  Figure  4  shows  how  I(x;e)  varies  with  x  for 

several  values  of  W,,  and  e  ■  0  and  e2  ■  0.5.  The  values  of  x  ,  x  ,  and 

5  m’  p* 

<x>,  and  the  corresponding  irradiances  1,1,  and  1  for  the  values  of  W, 

m  p  c  ) 

considered  are  noted  in  Table  2. 

The  variaace  of  the  aberration  p^cosQ  is  reduced  even  further  if  an 

3 

appropriate  amount  of  p  cos0  aberration  is  also  introduced.  For  a  given  value 
of  ,  the  appropriate  amounts  of  and  that  give  minimum  variance  may  be 
obtained  from  the  radial  Zernike  annular  polynomial  R*(p;e),  where7 

Rj(p;e)  * 

(56) 

Thus,  in  the  case  of  a  circular  pupil,  the  variance  is  reduced  by  a  factor  of 

3 

100  if  we  introduce  p  cos6  and  p  cos6  aberrations  with  ■  0.3W^  and 

W3  *  -1.2Wj.  Hence  the  peak  value  of  the  PSF  for  a  circular  pupil  aberrated 

by  a  small  value  of  W,  and  W,  ■  -1.2W-  occurs  at  x  ■  -0.6W.F.  According  to 

j  J  j  m  j 

Eq.  (43),  the  corresponding  centroid  occurs  at  <x>  ■  -0.4W3F.  Therefore,  the 


10(l-»-4e2>elt)p5  -  12(l+4e2+4e‘*-»e6)p3  +  3(l+4e2+i0el>+4e6+e8)f 

(l-e272  l(l+4e2+€l*)(l+9e2+9e,,+e6)  ]1/>2 
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Table  2.  Same  as  Table  1,  except  that  the  aberration  is  secondary  coma, 
W(h,0)  *  p5  cos0,  where  is  in  units  of  X. 


W5 

X 

m 

X 

P 

<x> 

I 

m 

I 

P 

I 

c 

1(0) 

0 

0 

0 

0 

1.0000 

1.0000 

1.0000 

1.0000 

(0) 

(0) 

(0) 

(1.0000) 

(1.0000) 

(1.0000) 

(1.0000) 

0.5 

0.50 

0.49 

1.00 

0.8150 

0.8153 

0.4114 

0.4955 

(0.63) 

(0.62) 

(1.75) 

(0.8400) 

(0.8402) 

(0.0760) 

(0.1768) 

1.0 

1.00 

0.83 

2.00 

0.4464 

0.4664 

0.0025 

0.2332 

(1.25) 

(1.21) 

(3.50) 

(0.4948) 

(0.4966) 

(0.0282) 

(0.0002) 

1.5 

1.50 

0.81 

3.00 

0.1685 

0.3237 

0.0098 

0.1873 

(1.88) 

(1.74) 

(5.25) 

(0.2003) 

(0.2196) 

(0.0130) 

(0.0009) 

2.0 

2.00 

1.11 

4.00 

0.0420 

0.2523 

0.0073 

0.1389 

(2.50) 

(1.71) 

(7.00) 

(0.0573) 

(0.1478) 

(0.0074) 

(0.0065) 

2.0 


separation  between  the  peak  and  the  centroid  is  0.2W^F.  For  large  values  of 
Wj,  minimization  of  variance  with  respect  to  and  does  not  lead  to 
a  maximum  of  the  PSF. 

As  an  example,  we  consider  the  PSF  aberrated  by  an  aberration 

W(h,0)  -  (W5P5  +  W3P3)  cose,  (57a 

where 

w3  *  -1.2  w5(l+4e2+4e*+e6)/(l+4e2+e%).  (57b 

According  to  Cq.  (56),  the  point  in  the  image  plane  with  respect  to  which  the 
aberration  variance  is  minimized  is  given  by 

x  »  -0.6  w,F(l+4e2+10e,*+4e6+e8)/(l+4e2+e,,).  (58) 

m  j 

Substituting  Eq.  (57)  into  Eq.  (19),  we  obtain  the  centroid 

<x>  -  -W5F(0.4+2E2♦7.2e,,+2e6+0.4e8)/(l+4e2+e,’).  (59) 

The  aberrated  PSF  along  the  x-axis  is  obtained  by  substituting  Eq.  (57)  into 
Eq.  (A6).  We  find  that  it  is  given  by  Eq.  (50),  where 

B  -  (2t2  W5  +  2tW3  -  xg)tl/2  .  (60) 

Figure  5  shows  the  aberrated  PSF  l(x;e)  for  several  values  of  with  W3  given 

by  Eq.  (57),  and  e  ■  0  and  e  “  0.5.  The  values  of  x  ,  x  and  <x>,  and  the  - 

m  p 

corresponding  irradiances  1,1,1  are  given  in  Table  3.  Note  that  x  ,  x  , 

tope  m  p 

and  <x>  are  all  negative.  Moreover,  their  magnitude  for  the  values  of 
considered,  especially  in  the  case  of  e2  “  0.5,  is  very  large.  Therefore,  in 
Figure  5,  the  horizontal  coordinate  is  chosen  to  be  x  -  x  . 


Table  3.  Sane  as  Table  1,  except  that  the  aberration  is  a  combination  of 
primary  and  secondary  coma  given  by  Eq.  (57). 


X 

X 

<x> 

I 

I 

I 

1(0) 

5 

m 

p 

m 

P 

c 

0 

0 

0 

0 

1.0000 

1.0000 

1.0000 

1.0000 

(0) 

(0) 

(  0) 

(1.0000) 

(1.0000) 

(1.0000) 

(1.0000) 

1.0 

-0.60 

-0.59 

-0.40 

0.9676 

0.9682 

0.8763 

0.3721 

(5.0) 

(-5.60) 

(-5.60) 

(-5.35) 

(0.8832) 

(0.8832) 

(0.0005) 

(0.0039) 

2.0 

-1.20 

-1.18 

-0.80 

0.8765 

0.8784 

0.5870 

0.0030 

(10.0) 

(-11.19) 

(-11.23) 

(-10.69) 

(0.6101) 

(0.6128) 

(0.0000) 

(0.0014) 

3,0 

-1.80 

-1.77 

-1.20 

0.7429 

0.7459 

0.2981 

0.0014 

(15.0) 

(-16.79) 

(-16.94) 

(-16.04) 

(0.3353) 

(0.3558) 

(0.0000) 

(0.0000) 

4.0 

-2.40 

-2.37 

-1.60 

0.5886 

0.5914 

0.1173 

0.0465 

(20.0) 

(-22.38) 

(-22.80) 

(-21.38) 

(0.1493) 

(0.2296) 

(0.0000) 

(0.0008) 
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Gaussian  Illumination 


As  an  example  of  a  radially  symmetric  illumination,  we  consider  a  Gaussian 
annular  pupil,  i.e.  one  for  which 

A(h)  -  Aq  exp[-y(h/a)2J 
-  Aq  exp(-yp2), 

where  Y  >  0.  The  aberration- free  PSF  for  such  a  pupil  may  be  obtained  by 
sunstituting  Eq.  (61)  into  Eq.  (1)  or  Eq.  (A6)  (Appendix  A]  and  letting 
W(h,0)  *  0.  Thus,  we  obtain 

l(r;Y;e)  «  (Y/(e”Y-  e~Y£  )]2  1(0, Y;e)  (/  exp(-Yt)  Jn(irr  t1^2)dt]2,  (62) 

C2  08 

where 

I(0;y;£)  -  (Tra2A^/XR)2  (e‘Y  -  e~Y£  )/y]2  (63) 

is  the  central  value  of  the  PSF. 

Primary  Coma 

If  we  let  $  ■  0  and  substitute  Eq.  (61)  into  Eq.  (A6),  we  find  that, 

along  the  x-axis,  the  aberrated  PSF  in  the  presence  of  primary  coma  of  the 
-  3 

type  W^p  cos0  may  be  written 


(61a) 

(61b) 


I(x;y;c)  *  (Y /(e”Y  -  e_Y£  )]2  1(0, y;e)  (/  exp(-Yt)  Jrt(irB)dt]2,  (64) 

e2  0 

where  B  is  given  by  Eq.  (51). 

Sunstituting  Eq.  (42)  with  n  *  3  and  Eq.  (61)  into  Eq.  (29),  we  find 
that  the  centroid  of  the  PSF  is  given  by 


(65) 


Prom  the  radial  polynomial  S^(p;e)  for  the  Gaussian  illumination7,  we  find 
that  the  point  in  the  image  plane  with  respect  to  which  the  aberration  variance 
across  the  Gaussian  pupil  is  minimized  is  given  by 


-2WF  2  •  ^  -  e~Y) 

m  3  Y  2 

e  (l+ye2)  ~  e~Y(l+y) 


For  small  values  of  W, ,  the  peak  value  of  the  PSF  occurs  at  x  . 

j  in 


We  now  consider  some  numerical  results  for  y  *  1,.  which  corresponds  to 

—2 

a  Gaussian  illumination  with  an  irradiance  of  e  at  the  edge  of  a  circular 

pupil  relative  to  the  irradiance  at  its  center.  Figure  6  shows  how  l(x;l;e) 

varies  with  x  for  several  values  of  and  e2  ■  0  and  0.5.  The  values  of 

x  ,  x  ,  and  <x>,  and  the  corresponding  irradiances  1,1,  and  1  for  these 
on  p  iq'  p '  c 

values  of  are  given  in  Table  4. 


Secondary  Coma 


The  aberrated  PSF  along  the  x-axis  in  the  presence  of  secondary  coma  of 
the  type  p5cos0  is  given  by  Eq.  (64)  where  B  is  given  by  Eq.  (55). 
Substituting  Eq.  (42)  with  n  *  5  and  Eq.  (61)  into  Eq.  (29),  we  find  that  the 
centroid  of  the  aberrated  PSF  is  given  by 


<„>  .  UW.P  «~2Te  1*1  i  Slh.  *  -«~2Y4_»  ILL*  UZtll  .  C67) 

e-2ye2_  e-2y 


Some  typical  numerical  results  are  presented  for  y  ■  1  in  Figure  7  and  Table  5. 


It  is  evident  from  the  data  given  in  Tables  4  and  5  that  the  centroids  of 
two  PSF' s  for  nonuniformly  illuminated  circular  pupils  aberrated  by  equal 
amounts  of  primary  coma  and  secondary  coma  are  different.  For. example,  when 
W^  *  Wj  *  IX,  <x>3  ”  1.37  and  <x>^  ”  1.12  where  <x>  is  in  units  of  XF. 

In  the  case  of  a  uniformly  illuminated  circular  pupil,  <x>^  *  <x>^  for 
as  may  be  seen  from  Tables  1  and  2.  Of  course,  for  an  annular  pupil,  the 
centroids  ^>3  «nd  for  W^  ■  are  also  different  whether  the  pupil  is 

uniformly  or  nonuniformly  illuminated. 
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Figure  6.  Seme  ea  Figure  1  except  that  the  amplitude  across  the  pupil  is 

Gaussian  given  by  Eq.  (61)  with  y  ■  l.  The  PSF's  are  normalised 
by  the  aberration-  free  value  l(0;l,e)  given  by  Eq.  (63). 


mwm 

Table  4.  Same  as  Table  1,  except  that  A(h)  is  a  Gaussian  given  by  Eq.  (61) 
with  y  *  1.  The  irradiances  given  here  are  normalized  by  the  aberration-free 
central  irradiance  l(0;l;e)  given  by  Eq.  (63). 


w3 

X 

m 

X 

P 

<x> 

I 

m 

1 

P 

I 

c 

1(0) 

0 

0 

0 

0 

1.0000 

1.0000 

1.0000 

1.0000 

(0) 

(0) 

(0) 

(1.0000) 

(1.0000) 

(1.0000) 

(1.0000) 

0.50 

0.61 

0.60 

0.69 

0.8805 

0.8806 

0.8670 

0.4567 

(0.76) 

(0.76) 

(1.42) 

(0.9288) 

(0.9288) 

(0.1126) 

(0.0602) 

1.00 

1.22 

1.15 

1.37 

0.6013 

0.6062 

0.5590 

0.1708 

(1.51) 

(1.51) 

(2.84) 

(0.7435) 

(0.7435) 

(0.1273) 

(0.0348) 

1.50 

1.82 

1.40 

2.06 

0.3205 

0.3672 

0.2479 

0.1199 

(2.27) 

(2.24) 

(4.25) 

(0.5112) 

(0.5122) 

(0.0014) 

(0.0033) 

2.00 

2.43 

1.46 

2.75 

0.1305 

0.2947 

0.0624 

0.0773 

(3.03) 

(2.93) 

(5.67) 

(0.3005) 

(0.3065) 

(0.0399) 

(0.0000) 
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V.  DISCUSSION  tiND  CONCLUSION 


We  have  defiaed  the  LOS  of  aa  optical  system  in  terms  of  the  centroid  of 
its  PSF.  Since  the  imaging  properties  of  an  optical  system  are  determined  by 
its  pupil  function,  the  centroid  of  the  PSF  is  no  exception.  By  expressing 
the  centroid  in  terms  of  the  pupil  function,  it  is  easy  to  show  that  the  wave 
diffraction  optics  and  ray  geometrical  optics  give  identical  expressions  for 
the  centroid,  regardless  of  the  shape  of  the  pupil  or  the  amplitude  and  phase 
distributions  across  it.  This  is  inspite  of  the  fact  that  the  two  PSF's  are 
quite  different  from  each  other.  Although  the  LOS  of  an  optical  system  can  be 
obtained  from  the  centroid  of  its  PSF,  or  from  the  slope  of  the  imaginary  part 
of  its  OTF  evaluated  at  the  origin,  in  optical  design  and  analysis,  the 
simplest  way  to  obtain  the  LOS  would  be  to  determine  the  centroid  of  the  ray 
spot  diagram.  The  idea  is  that,  since  ray  tracing  would  be  needed  to 
calculate  the  aberrations  of  the  system  any  way,  one  might  as  well  trace  the 
rays  up  to  the  image  plane  and  calculate  their  centroid  without  calculating 
the  diffraction  PSF.  Of  course,  the  precision  with  which  the  centroid  would 
be  calculated  would  depend  on  the  number  of  rays  used,  just  as  it  would  depend 
on  the  number  of  rays  used  to  calculate  the  aberration  function  and  the  number 
of  points  used  to  calculate  the  aberrated  PSF,  for  example,  using  an  FFT 
algorithm. 

In  the  case  of  an  aberrated  system  with  an  annular  pupil  and  a  radially 

symmetric  illumination,  e.g.,  Gaussian,  the  LOS  may  be  determined  from  the 

aberration  coefficients  of  the  appropriate  orthogonal  Zernike 

polynomials  S*(p;e)  cos0  and  S^(p:e)  sin6.  When  the  amplitude  across  the 
n  n  ^ 

pupil  is  uniform,  these  polynomials  reduce  to  polynomials  R  (p;c)  cos6  and 

1  n 
Rn(P;e)  sin6,  respectively.  The  results  for  a  circular  pupil  may  be  obtained 

by  letting  e  *  0.  In  practice,  however,  given  a  certain  optical  system,  the 

simplest  way  to  determine  the  centroid  would  be  to  use  its  point  spread 

function  as  measured,  for  example,  by  a  phctodetector  array. 

We  have  shown  that  for  a  uniformly  illuminated  pupil,  since  the  centroid 
of  an  aberrated  PSF  depends  only  on  the  aberration  along  the  perimeter  of  the 


-39- 


pupil,  aberrations  such  as  Wq  pacos6  and  Wq  R^(p;0)  cos6  (and  similarly 

W  p°sind  and  W  R^(p;0)  sind)  across  a  circular  pupil  with  different 
n  n  n 

aberration  orders  n,  give  aberrated  PSF's  with  centroids  that  depend  only  on 
the  value  of  Wq.  Thus,  for  a  given  value  of  Wq,  different  aberration  orders 
give  PSF's  with  the  same  centroid. 

We  have  obtained  numerical  results  on  the  PSF's  aberrated  by  primary  and 
secondary  coma  for  circular  as  well  as  annular  pupils  with  uniform  and 
Gaussian  illuminations.  In  the  case  of  uniform  illumination,  for  a  given 
value  of  W^,  the  peak  value,  is  higher  for  E2  -  0.5  than  for  E  ■  0.  The  peak 

value  as  well  as  the  centroid  occur  at  larger  value  of  x  when  £  is  no  zero 

compared  to  when  it  is  zero.  The  peak  value  corresponds  to  minimum  aberration 

variance  (I  ~  I  )  for  W_  <  1.5X  when  £  *  0  and  for  W_  <  2.5X  when  E2  ■  0.5. 

pm  j  ~  J  — 

Except  when  W_  ~  0,  I  is  quite  small  compared  to  I  .  As  W_  increases,  <x>  and 
j  c  p  J 

x  increase  linearly  with  it.  However,  x  first  increases  monotonically  with 
m  p 

Wj  and  then  fluctuates  in  a  series  of  maxima  and  minima.  For  £  *  0.5  the 

maxima  and  minima  are  widely  spaced  and  are  less  pronounced  compared  to  those 
for  £  ■  0. 


In  the  case  of  secondary  coma,  the  peak  value  does  not  occur  for  larger 
and  larger  values  of  x  as  increases,  i.e.,  xp  does  not  increase  monotonically 
with  Wj.  For  example,  when  £  ■  0,  *  0.83  for  ■  IX  and  x^  ■  0.81  for 

W,  ■  1.5X.  Similarly,  when  E2  ■  0.5,  x  ■  1.74  when  W-  ■  1.5X  and  x  *  1.71 


n  .  —  a  •  y  w  auia  ama  a  t  i  ^  v  •  ^  t  a  —  a  •  r  ~r  wuvu  n  p  a  •  a  at  >  a 

5  p  5  p 

when  W ^  *  2X.  As  in  the  case  of  primary  coma,  <x>  increases  with  £  for  a  given 


value  of  W^. 


In  the  case  of  Gaussian  illumination  numerical  results  are  obtained  for 


Y  *  l.  I  and  I  are  higher,  and  x  ,  x  and  <x>  are  smaller  for  Gaussian 
m  p  m  p 


illumination  than  those  for  uniform  (y  “  0)  illumination.  It  should  be  noted 


that  I  and  I  are  normalized  by  the  aberration-free  per.ic  irradiance,  and  it  is 
m  p 

only  these  normalized  values  Chat  are  higher  for  Gaussian  illumination.  [This 


normalization  is  different  for  Gaussian  and  uniform  illuminations  as  may  be 
seen  by  comparing  Eqs.  (41)  and  (63)].  Otherwise,  for  a  given  circular  or 
annular  aperture  and  for  a  given  total  power  E,  an  unaberrated  pupil  gives 
maximum  central  irradiance  when  it  is  uniformly  illuminated. 12 


The  results  given  here  are  applicable  to  both  imaging  systems,  e.g., 
those  used  for  optical  surveillance,  as  well  as  to  laser  transmitters  used  for 
active  illumination  of  a  target.  In  both  cases,  the  LOS  of  the  optical  system 
is  extremely  important.  A  LOS  error  of  a  surveillance  system  will  produce  an 
error  in  the  location  of  the  target.  In  the  case  of  a  laser  transmitter,  a 
large  LOS  error  may  cause  the  laser  beam  to  miss  the  target  altogether. 

Whereas  for  static  aberrations,  we  may  be  able  to  calibrate  the  LOS,  for 
dynamic  aberrations-  it  is  the  analysis  given  here  that  will  determine  toler¬ 
ances  on  aberrations  of  the  type  pncos6  and  pnsin6  in  the  case  of  annular  or 
circular  pupils. 

% 

Although  we  have  defined  the  LOS  of  an  optical  system  in  terms  of  the 
centroid  of  its  PSF,  it  could  have  been  defined  in  terms  of  the  peak  of  the 
PSF  (assuming  that  the  aberrations  are  small  enough  so  that  the  PSF  has  a 
distinguishable  peak).  As  pointed  out  in  Section  II,  for  an  aberration-free 
PSF,  its  peak  value  and  its  centroid  both  lie  at  its  origin,  regardless  of  the 
amplitude  variations  across  its  pupil.  The  two  are  not  coincident  when 
aberrations  are  present.  The  precise  definition  of  the  LOS  will  perhaps 
depend  on  the  nature  of  the  application  of  the  optical  system.'  Moreover,  in 
practice,  only  a  finite  central  portion  of  the  PSF  will  be  sampled  to  measure 
its  centroid,  and  the  precision  of  this  measurement  will  be  limited  by  the 
noise  characteristics  of  the  photodetector  array. 
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APPENDIX.  PSF  for  an  Annular  Pupil  with  Radially  Symmetric  Illumination  and 
Coma  Aberration 

Consider  an  optical  system  having  an  annular  pupil  with  radially 
symmetric  illumination  A(h)  and  coma  aberration 

W(h,6)  ■  l  pQ(W  cos0  +  W'  sin6). 

n  n 

n 

where,  as  in  Eq.  (25),  nhe  prime  on  the  summation  sign  indicates 
over  odd  integral  values  of  n.  Note  that  h  *  pa  and  e  £  p  l. 
coordinates,  Eq.  (1)  for  the  aberrated  PSF  may  be  written 

l  2  it 

I(r,*;e)  -  (a2/XR)2|/  /  A(h)  exp{Tril2W(h,9)  -  r  p  cos(0-*)J}  p  dp  d0|2,  (A2) 

e  0  8 

where  we  have  assumed  that  W(h,6)  and,  therefore,  W  and  W'  are  in  units 

n  n 

of  X.  Integration  over  0  in  Eq.  (A2)  can  be  carried  out  if  we  let 


(Al) 

a  summation 
In  polar 


2W(h,0)  -  r^p  cos( 0-*)  ■  B  cos(0-*),  (A3) 

where 

B2  -  ( l '  2  Wq  pn  -  r#  p  cos*)2  +  (['  2  W'  p"  -  r#  p  sin*)2  (A4) 

n  n 

and 

tan*  -  (l'2W'pn  -  r  sin*)/(^2W  p“  -  r  cos*)  (A5) 

q  8  n  8 

n  a 

Thus,  Eq.  (A2)  becomes 

1 

I(r,*;e)  -  (7Ta2/XR)2  {/  A(h)  J.(TfB)dt]2  ,  (A6) 

e2  0 

where  we  have  let 


(A7) 


?.L> 
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